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ABSTRACT 


A Fortran IV computer program provides convenient simulation of an all- 
digital phase-lock loop (DPLL). The DP LL functions as the primary phase processor for 
the Ohio University Omega prototype navigation receiver. 

This project has been supported by NASA Langley Research Center Grant 
NGR 36-009-017 for application to general aviation VLF receiver systems. 

II INTRODUCTION 

This paper describes a Fortran IV simulation study of the all-digital phase— lock 
loop (DPLL) which was designed by the NASA Omega staff of the Ohio University Avionics 
Engineering Center. The DPLL forms the heart of the Omega navigation receiver prototype. 
Through the DPLL, the phase of the 10.2 KHz Omega signal is estimated when the true 
signal phase is contaminated with noise. The study has provided a convenient means of 
evaluating loop performance in a variety of noise environments, and has proved to be a 
useful tool for evaluating design changes. 

The goals of the simulation were threefold: 

A. To analyze the circuit on a bit-by-bit level in order to evaluate the overall 
design; 

B. To see easily the effects of proposed design changes prior to actual breadboard- 
ing; and 

C. To determine the optimum integration time for the DPLL in an environment 
typical of general aviation conditions. 

III DESCRIPTION OF THE SIMULATION 

The digital nature of the circuit indicated the use of the digital computer as 
the most convenient simulation tool. Since the circuit was already built and was, in fact, 
yielding phase information which compared favorably with the Tracor series 599R receiver 
used for control-group data collection purposes, the circuit model was wholly deterministic. 

The generation of zero-crossing information for the Omega signal in noise, 
however, was necessarily stochastic tn nature. A complete description of the signal- 
generating procedures is given in part VI. 

CORDET is based on the technique of periodic scanning; that is, each portion 
of the DPLL is observed and updated with each clock pulse, significant information about 
the state of the system is recorded, and the process is repeated. In the case of CORDET, the 
significant information includes the states of all logic circuit elements, including phase 
difference information referenced to a local oscillator. The phase information is,of course, 
of primary importance in an Omega navigation receiver. The ability to compare the DPLL 



output with the true phase of the simulated Omega signal/ and to provide graphical output, 
are advantages of a computer simulation over a laboratory hardware experiment. A user's 
technical description is provided in the Appendix. 

IV DESCRIPTION OF CIRCUIT 

The DPLL follows closely the ideas of J.M. Clark^ ^ , first presented in 1968. 
The DP LL presently in use at Ohio University is essentially a highly adapted version of his 
digital phase tracking filter. A detailed discussion of the loop circuit operation is given 

by Chamberlin^. 

The circuit employs a 6-bit counter which cycles once every 64 pulses of the 
652.8 KHz clock (see Figure 1). This counter therefore cycles at a rate of once each 10.2 
KHz (=652.8 KHz -*■ 64). The number in the counter is compared continuously to the six 
most significant bits in the up-down counter, which Is in turn fed by the phase detector. 
During the brief period when the numbers in the two counters are bit-wise identical, the out 
put of the comparator goes to the logical 1 state. Although the duration of this pulse 
may vary, the uni-directional counter generally advances at a much faster rate, indicating 
a pulse width of 1/652.8 KHz, or approximately 15 psec. This pulse is used to trigger a 
monostable, creating a window whenever the contents of the two registers are identical. 

When viewed through this window, the zero crossings of the incoming signal are 
seen to be either leading or lagging in relative phase. When a lag is indicated, the up- 
down counter receives a down-count command, and the next window is created earlier In 
time. In this way the edge of the window gating function is brought more and more closely 
into coincidence with the zero crossings of the true Omega signal. Similarly, when a phase 
lead is indicated, count-up commands cause the window to be created later in time. 

V THE DPLL AS A CROSS-CORRELATION DEVICE, . 

In a cross-correlation type receiver scheme, two signals are combined to give 
a d.c. level by first delaying one signal In time, multiplying their resulting, signals, and 
integrating or lowpass filtering (see Figure 2). To see how the DPLL performs this same 
function in an all-digital sense, it is first necessary to view the circuit operation under no- 
noise conditions. Under these conditions there is no phase shift due to noise perturbations, 
nor is there a Doppler shift. The only phase error is then due to initial phase difference at 
turn-on. 

Assume that the initial phase difference is such that a count-up command is 
given to the up-down counter. If N bits are being integrated, there will be a delay of 2 
up-count commands before the comparator logic "sees" a difference in time between suc- 
cessive occurrences of identical numbers in the two counters. The time delay in such a 
case is (2 n /10.2 x 10^) seconds. 

The output of the comparator, 8, is the estimate of the Omega signal phase, 

0q , When lock is obtained, 9 = 9 q, and the binary number contained in the six most sig- 
nificant bits of the up-down counter Is the phase difference, analogous to the d.c. level 
obtained In an analog cross-correlation device. 
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VI SIGNAL INPUT SIMULATION 


Since the receiver front end provides a 15 Hz bandwidth narraw-band filter, 
the derived phase information can be considered to result from a noise-plus-signal input 
to a narrow-band system. 

The composite input signal can be expressed as 

g (t) — f (t) + A cos (<*» t + t) (1) 

where u m is the center frequency of the filter, and f (t) Is the narrow-band noise process 

f(t) = E (t) co5:(u m t + 0(t)). (2) 


[3] 

It can be shown ■ that the amplitude process E(t) is Rayleigh distributed, while the phase 
process 0(t) is uniformly distributed on the interval (0, 2tt). 

Now, f(t) can be written in the form defining its quadrature components, 

f(t) = e (t)cosw t + e (t) sin w t (3) 

c ms m 

[ 31 ' 

where e (t) and e (t) are independent Gaussian processes 
From equation 1, g (t) may be expressed as 


g (t) - f(t) + A cos tycosu t + A sin t sin w t. (4) 

m m 


Substituting for f(t), and using 

cos (a + b) — cos a cos b - sin a sin b. 


we have 

g (t) = (e (t) + A cost) cosu t + (e (t) + A sin t) sinu t. (5) 
c m s m 

Equation 5 may be put in the form 

g(t) « R(t) cosfw^t + 9 (t) ) (6) 

where R(t) = { [e (t) + Acos^l^ + [e (t) + Asintl^} 5 (7) 

c s 
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and 


e (t) + A sin 

9 (t) = arctan — L - , ..-..- — — (8) 

w e (t) + A cos $ 

c 

Equation 8 is therefore the input phase representation to the phase detector, 

where e (t) and e (t) are statistically-independent Gaussian white noise processes with 
c s 

zero mean and unit variance, and A is the normalized signal amplitude. In general, 
can be a function of time to represent the Doppler shift in the dynamic case. At speeds 
less than Mach I, however, the Doppler shift is well" below the 5.06° (1/64 of a 10.2 KHz 
cycle) quantizing interval of the DPLL, and so ijican be considered constant over the 
sampling interval. 

A block diagram of the noise input phase simulation is given in Figure 3. 

VII CONCLUSIONS 

In determining the optimal integration time of the first order loop, consider the 
case of a small aircraft in a noisy environment. As the time-multiplexed Omega transmission 
sequence progresses, the aircraft will change position between successive transmitting time 
intervals of the 10.2 KHz frequency. For an aircraft flying a course coincident with a 
station pair baseline (worst case), the change in position will amount to 12° of 10.2 KHz 
phase at a speed of 200 mph, conservative for private aircraft. For this reason, step phase 
inputs of 12° were taken both in a noise environment of 10 decibels SNR and in the no-noise 
case. 


To determine optimal integration time, phase error at the end of the 625 msec, 
sampling interval is plotted vs. integration time in Figures 4 and 5. As could be expected, 
the best Integration time in noise is one which is as long as possible without leaving a 
residual phase error at the end of a sampling interval in the no-noise case. For aircraft at 
200 mph, 0.1 sec., or 10 bits of integration proves to be optimal. 

VIII SUMMARY 

The Fortran program CORDET can prove to be a useful tool in analysis and design 
of a first order all-digital phase-lock loop. The DPLL simulated in CORDET is suggested as 
useful for designers of digital processing Omega receivers or other VLF signal processing 
systems. 

IX ACKNOWLEDGEMENTS 

Research on Omega at Ohio University is funded by NASA. This work is being 
done at the Avionics Engineering Center headed by Dr. Richard H. McFarland. Mr. Ralph 
Burhans is project engineer and Dr. R. W. Lilley is staff consultant on the Omega project. 
The author wishes to thank these persons for their assistance in the preparation of this paper. 
Additional assistance was provided by K. A. Chamberlin and J. E. Grover, graduate 







Phase Error, Phase Error, 

Degrees of Lag Degrees of Lag 



Bits of Integration 

Figure 4. Phase Error at End of Sampling Interval vs. Bits of 
Integration, No-Noise Case. 
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Figure 5. Typical Plot of Phase Error at End of Sampling Interval 
vs. Bits of Integration, 10 Decibel SNR. 
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APPENDIX: USER'S TECHNICAL DESCRIPTION 
OF FORTRAN PROGRAM CORDET 


A. Genera! Description, 


CORDET is a logic simulation of the all-digital phase-lock loop (DPLL) presently 
employed in the Omega receiver prototype. 

Output provides graphic display of the timing diagrams of the circuit, plots of 
phase difference vs. time, and a histogram of phase difference between the DPLL output and 
the object-lock (input) frequency. 

Provisions have been made to vary the number of bits of integration, the signal- 
-to-noise ratio of the input signal, and-the initial phase difference between the DPLL output 
and the object- lock frequency. 

B. Capabilities and Limitations . 

1. The program assumes the presence of a local oscillator in the circuit with 
a frequency of 64 times the object- lock frequency. 

2. Compile and load time is 15 seconds. Execute time is approximately 
0.15 seconds per cycle of object-lock frequency specified (when only phase difference and 
histogram are plotted). 

3 

3. - Approximately 21 x 10 bytes of storage are used, 

4. Subroutines HISTG and RANDNR are called from the Ohio University 
computer center request-call library. RANDNR uses the multiplicative congruential method 
of pseudo-random number generation. HISTG is a plotting subroutine. 


Technical Description. 


Language: 

FORTRAN IV 

Number of arguments: 

11 

1 . Data card 1: 

NCYC, MNUMB, NNUMB, DCYC, 
DELPHI, NOISE, A 

Mode: 

INTEGER (1, 2, 3, 6), REAL (4, 5, 7) 

2. Data card 2: 

TD, PHA, HISTOG 

Mode: 

INTEGER 

3. Data card 3: 

TITLE 

Mode: 

ALPHANUMERIC 
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Argument definitions: 


NCYC - is the number of object-frequency cycles to be 
considered. 


MNUMB 

NNUMB 

DCYC 

DELPHI 

NOISE 

A 

TD 

PHA 

HISTOG 

TITLE 


- is the number of bits in the up-counter. 

- is the total number of bits in the up-down 
counter. 

- is the duty cycle of the monostable. 

- is the initial phase difference in degrees of 
lag (DPLL output lags input). 

- is zero for no-noise case, 1 otherwise. 

- is voltage signal-to-noise ratio in rational 
number form (not db). 

- is 1 if timing diagrams are desired, 0 otherwise. 

- is 1 if phase difference vs. time plots are 
desired, 0 otherwise. 

- is J_ if histogram is desired, 0 otherwise. 

- is title for histogram plot. 


Figure A-l represents the positive- logic implementation of the DPLL. CORDET 
performs one pass through a DO- loop for each pulse of the 64 x 10.2 KHz clock. New values 
for the logical variables indicated in the diagram are assigned on each pass. 

Changes in the DPLL circuit will of course necessitate changes in the logical 
statements in CORDET d s main program. These changes are facilitated by liberal use of com- 
ment cards in the program deck. 



LMONO 


Figure A-l. Positive-Logic Implementation of the DPLL. 
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c INITIATE AND DIMENSION 

INTEGERS REF, S( 24) , INDI *NS, TO ,PHA ,HTSTOG 
INTEGER PHADIF,PH(650) 

LOG I CAL LMONO.UP, I NDLCJC * UPPE * ON , L OM EG A • I ND, LMONO 2 * EX * POMEGA- 
C -OM ' AND ' AM "MUST BE' D I MENS I CM ED TO 2* ( 2* ^MNUMB ) +1 1 


LOGICAL^ 1 OM{ 129) , AM ( 1 29 ) ' 

DIMENSION ARRAY (650) 

EX = .FALSE* 

LM0N02 = .FALSE, ' ' ‘ " . 

INDLOC = .FALSE. 

POMEGA = .FALSE. 

READ (1,133 NCYC • MNUMB, NNUM3 , DCYC , DELPHI , NO l SE, A 


13 ‘ FORMAT (3.15, 2F10.2*T2,F5.2) 7 

READ (1,14) TO,PHA ,HISTOG 

14 FORMAT (315) * 

C 7 THETA IS INITIAL PHASE SH I FT • -2*Pi *DELPH I /360 * 


• THETA - DELPHI *0 .0 174532 9 
SINE = SIN('THETA) . 

COSINE = COS ( THETA) 

DO 57 1=1,129 

' ~OM< I ) ~i F ALSE. - '” 

AM(.I )= .FALSE. . 

57 CONTINUE 

JJ = 0 

JJJ ■= 0 

11 = 0 . 

MONO = 0 
M = 0 

— • N = 0 • “■ - ' ‘ 

MR EG M — 2**MNUMB- 
MREGN = 2**NNUM9 
HALF = FLOAT (MREGM 3/2 . 

MHALF =- HALF ' “ 

K = MREGM 

NCYC P •= NCYC^MREGM 

NCYCT = NC Y C * * 1 

X = ~I DELPHI /3 60 .') *MREGM ’ “ 

IX = X 

III = 0 

ji r • = 64 

TNDCH *=• " 0" • 


C DO-LOOP BEGINS WITH OCCURRENCE OF FIRST.. CLOCK PULSE 

DO 12 1=1, NCYC P * 

C M IS NUMBER STORED IN UP-CGUNTER — " . 

IF (M-MREGM) 10*11,11 

1 1 M = 0 * - 

10 CONTINUE ' * 

C ; B INARY 'NUMBER IN UP— ON* COUNTER’ IS "TRUNCATED DUE"' TO BIT''" 

C- —INTEGRATION. NTRUN IS RETURNED . 

CALL TRUN ( N NUMB ,MNUM3 , N ♦ N TR UN , 3 3 

c IF m- NT RUN, REF GOES HIGH 

"IF ( M-NTRUN) 1,2,1" L - 

1 ■ REF = 0 . . 

GO TO 3 •' 

2 REF = 1 ‘ - 

3 CONTINUE - ■ -- — 

C MONOS PROVIDES OUTPUT OF MONOSTABLF. K IS TIMER FOR MONO. 

C LMONO IS LOGICAL OUTPUT OF MONOSTABLE. 

K = K + 1 

• ' CALL" MONOS "( K , REF, DCYC, MREGM , LMONO ) . r — - - 

C : — L OMEGA IS OMEGA SIGNAL ZERO GROSSINGS. 

c A IS. VOLTAGE SNR. POMEGA IS TRUE PHASE OF OMEGA SIGNAL. 

IF (IX.EQ.M) GO TO 30 ’ 

POMEGA =. FALSE. : ’ ■“ ‘ : " • 

GO TO 31 

30 POMEGA = .TRUE. 

31 CONTINUE 

IF ( NO I SE. EO .0 ) GO TO‘32 - 

IF (POMEGA) GO TO 34 • -■ 

GO TO 51 

. 34 I NDCM - — 1 # I NDCH+ 1 . 

IF ( INDCH. ECU1 .) GO TO 37 ‘ ‘ — ■ r - 

DO 3 8 UK = 1*129 . 

AM { I JK) =. FALSE. 

30 CONTINUE 

• C NOISE IS ADDED TO QUADRATURE COMPONENTS' OF : OMEGA SI GNAL • 

39 CALL ANORM (AN) 

ANUM = A*SINE+AN 

■CALL ANORM (AN) - 16 - .. ' •. . - . 



50 


37 

52 


51 


GO ‘TO '50 
=«TRUE» 

= .TRUE. 


DENOM - A*COSlNE+AN 
THETA I = AT AN2 . ( ANUM, DENCM > 
PCYC •= ( THETA I /6 .283 1 3) *MREGM 
INDEX = PCYC 
IE" { INDCH. EO. 1 )' 

AM ' ( MREGM+ INDEX) 

GO TO 51 

OM { MREGM+ INDEX) 

GO TO 51 
CO 5 2 I JK-1 ,129 
QM(IJK) =• FALSE. 

CONTINUE 
‘GO TO 3 9 
CONTINUE 
J I I - J I I + 1 
III = I I I TT 
IE ( I I I «EQ. 1 29 
IF( J I I *EQ. 129 ) 

LOMEGA 
IF ( ( I 


III = I. 

II = 1 

“ GM.( I I I ) .OR, AM { J I I ) 

EO .M) . AND . < POM EG A) ) L OMEGA - 


POMEGA 


32" 
C 


C 


c 

8 

9 

C 


I F ' ( NO I5E .EQ.O ) LOMEGA = POM EG A ’ • 

LMONO AND IOMEGA ENTER PHASE DETECTOR 

UP - LMO NO. A ND.L OMEGA 

ON — < • NOT * LMONO) , AND.LOMEGA 

• — 'PHASE - DETECT OR OUTPUT FEEOS'UP-DN COUNTER r 

.IF < UP ) N = N+ 1 
IF ( DN } N = N-l 

COUNT REGISTER MUST CONTAIN POSITIVE NUMBER 

IF< Ni> 3*9,9 ■“ 

N s MREGN+N 
CONT INUE 

UP-DN- COUNTER COUNTS MODULO-MREGN 



•*" IF { MREGN-N ) 23, 24 * 

24 

23 

N - 0 


24 

CONTINUE 

IF (TO.EQ.O) GO TO 
™ TF~ ( I— 1 ) 4 0,40,70*” 

50 1 



40 

CONTINUE 

TIMING DIAGRAM 

I S 


PLOTTED 


50 O' 


70 


50 1 


( ITN;NTRUN,M9 ref, lomegaylmono ,lmono2 vdn, up; indlocvex 


GO TO 12 


23 


22 


21 

20 

45 

c — 


65 

12 


502 
C 


27 


146 


PRINT 500 

‘ FORM AT ( 1 1 *Y2X , * I • , 3X * N* ,3X , » NTRUN* 
* , • POMEGA* ,4X , * UP • , 6X , • DN * , 4X , * LMONO 
*X»//) 

CONT INUE 
-"CALL PLOT! 

* M EGA) ’ ' 

CONT INUE 

IF ( ;(PHA.EO.O),AND.(HISTOG.EG.O) ) 

IF" ( .NOT. POMEGA) GO TO 65 
JJ ~ JJ + 1 

IF < JJ .EQ. 1 ) GO TO 2 8 
A J = JJ/10* 

■ i j = aj • ‘ ' — “■ ‘ • 

IF ( 1J.NE.AJ) GO TO 65 
.CONTINUE 

IF (M.LT.MHALFJ GO TO 20 

’ I F" ( NT RUN .LT .MHALF ) GO TO: 21 

P HA OIF - 2*{NTRUN-M) 

GO TO 45 

PH AD I F - < NTRUN +• MREGM-M)«2 

- GO TO 4 5 * ■ * 

IF ( NTRUN.LT. MHALF) GO TO 22 
PH AD I F = NTRUN— MPEG M— Mi) 

■CONT INUE 


‘3 X . 

, 3 X . 


M * , 4X, •REF* ,3X.‘* LOMEGA-* ■ 
» LMONQ2 * , 2X » * INDLOC ’ ,AX 


TO 50 2 


— ARRAY’ IS" 

JJJ = JJJ+1 
PH(JJJ) = PH AD IF 
CONTINUE 

-CONTINUE 

IF (PHA.EQ.O) GO 
CALL PLOT (PH, JJJ) 

CONT INUE 

" “IF ' ( HrSTQGVEQi'O ) 'GO 

ARRAY IS LOADED 

DO 27 I = 1 t JJJ 

ARRAY (ID = FLOAT ( PH ( I)) 

" CONTINUE 

DIMENSION TITLE! 20 ) 

READ (1,146) T I TLE 
FORM AT ( 20 A-4 ) - 


LOADED’ FOR" PLOT "SUBROUT IN E- — 


TO 503 
FOR HISTG 


SUBROUTINE 


- 17 - 


2X’ 

•E 

PO“ 



50 3 


C 

C 


9 


i 


2 


3 

4 

C 


_ rr 
c - 


1 o 


1 o 

• i i 

12 
1 9 


1 4 
1 8 

15 

1 6 
l 7 


OATA AST/* £ * / 

HISTOGRAM- IS PLOTTED 

CALL HISTG (JJJ t ARRAY* TI TL E . 1 » A ST ) 

CONT INUE 

" STOP ■ " 

END ... ■ • 

SUBROUTINE TRUN { NNU MB , MNU M3 * N , NT RUN * 8 ) ' 

THIS SUBROUTINE CONVERTS DECIMAL INTEGERS TO BINARY FORM 
AND TRUNCATES, SHIFTS D OWN " REMA I N ING BITS AND RETURNS A 
INTEGERS 8 < NNUMS ) . 

DO 9 I = 1 ► NNUMB 
B ( I ) = 0 

CONTINUE ’. T — - ’ 

1=0 

DUM = N 
DUM=DUM/2. 

I DUM-DUM 

IF (DUM- I DUM) 2,3*2 
EH I > = 1 

' DUM = DUM- V 5* ■ - — ~ — 

I F ( DUM ) 4, 4, 1 

B ( I ) = 0 

IF ( DUM ) 4,4,1 

’CONTINUE - - - — - 

THE BINARY NUMBER IS TRUNCATED 

DO 1 1 I - 1 * MNUMB 

B{ I ) = B( I + (NNUNB-MNUM8).J 

CONTINUE ' — "• •— ■ - - 

THE TRUNCATED BINARY NUMBER IS CONVERTED BACK TO DECIMAL 

NT RUN ;= 0 

DO 10 T — 1 , MNU MB 

K=Bcr ) *< 2 **-{ r-T ) i • — •' ' — - - 

NT RUN ■= NT RUN + K 

CONTINUE 

RETURN 


SUBROUTINE MONOS ( K , REF , DC YC, MRFGM , LMONO } r 
LOGICAL 'LMONO 
INTEGER^2 REF 

X = DCYC*MREGNI - 1 . 1 ' " - 

IX - X 

IF ( REF ) 10,10*12 

I F ( K- I X ) 1 4 • I 4 , I 1 

MONO = 0 *:■ - - - 

GO TO 18 1 ’ 

IF (K-IX) 14,-1 4, 19 • 

MONO = 1 

' K=0 ' - ; - ’ ■ 

GO TO 18 
MONO = 1 

IF ( MONO > 1 5 , 1 5 * 1 6 

"LMONO = . FALSE. * ’ - - — -- 

GO TO 17 
LMONO = « TRUE • 

RETURN 

END. - V • 


DEC I 


- 18 - 



SUBROUT I ME PLOT (PH,JJJ) - 

REAL LINEU29) 

INTEGER PHADIF,PH(JJJ) 

DATA DOT/* *»/, PLUS/ 9 + */, BLK/ * *./ , VL I NE/» j *V , AST/ * * • / ' ' 

DO 1 2 J J = 1 , J J J 
IF (JJ.ME.l) GO TO 5 
PRINT 3 

8 FORMAT ( * 1 * , 4 5 X , 'PHASE DIFFERENCE IN DEGREES VS, TIME IN MSEC# 

4 PRINT 2 

2 FORM AT ( 3X , » - l 80 * , 2 QX , » - 9 0 30X , * 0 * » 31 X , * 90 » , 28 X , * 1 80 * > 

6 PRINT 7 

7 FORMAT ( 4 X * * + * , 1 6 ( ’ + *)) 

5. CONTINUE 

NSP. = PHCJJ) + 65 
IF { NSP * G T #129) GO TO 15 

IF. < NSP.LT, 1 ) GO TO 15 

GO TO 3 

15 CONTINUE 

PRINT 2 3 , P H ( J J > 

2.3 FORMAT { 1 OOX * • PH= - * ,110) 

RETURN 

3. CONTINUE 

C- CLEAR LINE 



DO 1 J = 1,129 

UNE(J) - BLK 
1. CONTINUE 

L I NE < 1 ) as VL INE 

L INE (12 9) VL INE 

A I = J J / 1 0 * 

I I - A I 

I F { A I • NE . I I ) GO TO 9 
L I NE <17} ■ = DOT 
LI NEC 33) ~ PLUS 

L I NE { 49) = DOT 

LINE (65) - PLUS 

’ ’ ' LINE ( 81 ) — 'DOT 

LINE (97) - PLUS 

LINE { l 13 ) = DOT ‘ 


9 CONTINUE 


10 

‘LI NEC NSP) = AST 
IF ( A I * EQ « I I ) GO 
WR I TE< 3, 1 1 > LI NE 

TO 2 0 

1 1 

FORM AT { 4X , l 29 A 1 ) 


— 

GO TO 1 2 


20 

TIME = ( FLO A T( J J) * 
WR I TE (3,21) LINE 

# CO 9 30 392 16) 

21 

FORMAT ( 4X , 1 29 A 1 ) 
WRITE (3,22) TIME 


22 

FORMAT (*+*»4X»E8 

#2 ) 

1 2 

CONTINUE 

RETURN 



END 
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SUB ROUT I ME °LOTl { I * N , NT RUN ♦ M *REF • LOM E'GA. i_MO.NO ,LM0N02 * ON , UP * INOLOC * 
*EX .POMEGA) 

INTEGER* 2 REF 

LOGICAL ' LMONO. UP* I NDLQC.DN.LOMEGA. LM0NQ2 * EX. POMEGA 


INTEGER’ LINE(56) 

INTEGER BLK/» » / , VL T NE/ * \ * / 

DO l J = 1*66 

LINE{ J)"- BLK 

1 .CONTINUE 

LI NEC I ) = VL INE 

.00 2 K - 1 ,8 

LINE (8 *K> = VL I NE ' 

2 CONTINUE 

IE CREE.NE.l) GO. TD 3 
LINEC3) “ VLINE 
‘ " LINE ( l ) - BLK 

3 IF { *NOT *L OMEGA) GO TO 4 
LINE CIO) = VLINE 

LINE C 8 ) = BLK 

-4 “ IF C.NOT.POMEGA) GO TO 5 

L I NE ( l 8 ) = VLINE 
LINE { 16)- BLK 

5 IF(.NQT.UP) GO TO' 6 

LI NE (26)— VLINE ‘‘ ■ 

LINE (24)= BLK 

6 IF(.NOT.ON) GO TQ 7. 

LINE(34)= VLINE 

LINE ( 32 )= BLK' ' ‘ 

7 IF (.NOT.LMONO) GO TO 8 

L I NE ( 42 ) = VL INE 

LINE (40)= BLK 

~"~3 "-"IF ■r.NOT.LMONO'2) ‘GO“TO ' 9 ~ ‘ * 

LI NEC 50)= VLINE 
LINE ( 48 )= . BLK 

9 IF { »NQT • I NDLOC ) GO TO 10 

• L 1 NE'C'5 8 ) VL TN E * 

LINE ( 56 )= BLK 

10 IF (.NOT. EX) .GO TO 11 
LINE(66)= VLINE 

.. ... ' LINE" (64)= BLK “ ‘ 

11 WRITE(3.12) l . N* NTRUN.M.LI NE 

12 FORMAT { 1 X.2 I 4 . I 6.* 15 .5X.66A1 ) 
RETURN 

. END - • ’ — • 

100 -6.8 0.50 0.01 

Oil 

PHASE DIFFERENCE IN 64THS OF A CVCLE 
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